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ABSTRACT 

The Adaptive Optics Facility (AOF) project envisages transforming one of the VLT 
units into an adaptive telescope and providing its ESO (European Southern Observa¬ 
tory) second generation instruments with turbulence corrected wavefronts. For MUSE 
and HAWK-I this correction will be achieved through the GALACSI and GRAAL AO 
modules working in conjunction with a 1170 actuators Deformable Secondary Mirror 
(DSM) and the new Laser Guide Star Facility (4LGSF). Multiple wavefront sensors 
will enable GLAO and LTAO capabilities, whose performance can greatly benefit from 
a knowledge about the stratification of the turbulence in the atmosphere. This work, 
totally based on end-to-end simulations, describes the validation tests conducted on 
a C"^ profiler adapted for the AOF specihcations. Because an absolute profile calibra¬ 
tion is strongly dependent on a reliable knowledge of turbulence parameters tq and 
Lq, the tests presented here refer only to normalized output prohles. Uncertainties in 
the input parameters inherent to the code are tested as well as the profiler response 
to different turbulence distributions. It adopts a correction for the unseen turbulence, 
critical for the GRAAL mode, and highlights the effects of masking out parts of the 
corrected wavefront on the results. Simulations of data with typical turbulence profiles 
from Paranal were input to the prohler, showing that it is possible to identify reliably 
the input features for all the AOF modes. 

Key words: instrumentation: adaptive optics - methods: data analysis - atmospheric 
effects - site testing 


1 INTRODUCTION 

A reliable knowledge of the atmospheric turbulence distribu¬ 
tion allows for the setting of constraints on the design of new 
adaptive optics (AO) systems, for optimizing their perfor¬ 
mance and for making a prior assessment of their correction 
limitations. Such statements are true for AO modes such 
as laser tomography or multiobject adaptive optics (LTAO 
and MOAO, respectively) which require an optimization of 
the wavefront reconstruction along one or more determined 
lines of sight. In single-, multiconjugated and ground layer 
adaptive optics (SCAO, MCAO and GLAO, respectively), 
the knowledge of the turbulence distribution enables one to 
determine the level of correction expected for the system, 
its anisoplanatic behaviour and the residual error that will 
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impact as a degradation of the science point-spread function 
(PSF) in the telescope focal plane. 

Many techniques exist for acquiring information on 
the turbulence, such as balloon borne experiments to mea¬ 
sure temperature structure coefficients, SCIDAR, SLODAR, 
MASS-DIMM, etc (Lombardi, Navarrete & Sarazin 2014). 
Monitoring programmes have been carried out in Paranal 
for several years to increase the database that will enable 
to set a programme to optimize operations - see e.g. Mas- 
ciadri, Lascaux & Fini (2013) - and to provide parameter 
constraints for the next generation of VLT and even E-ELT 
instrumentation (Lombardi et al. 2013; Sarazin et al. 2013). 

The Adaptive Optics Facility (AOF) is an AO-oriented 
upgrade to be commissioned in one of the telescope units 
of the Paranal Observatory in 2016. It counts on the largest 
Deformable Secondary Mirror (DSM) ever built (1170 voice- 
coil actuators distributed behind a thin shell of 1120 mm of 
diameter) and on four Sodium laser side-launch telescopes 
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which are part of the new Laser Guide Star Facility. The 
AOF will operate in 3 distinct modes (SCAO, GLAO & 
LTAO) in accordance to the instruments attached to the 
2 Nasmyth ports (GRAAL+HAWK-I, GALACSI+MUSE) 
and to the Cassegrain port (ERIS) of the telescope. Each 
of the four Shack-Hartmann wavefront sensors counts on a 
40x40 squared subaperture (SA) array, out of which 1240 
SAs are considered valid and 1152 are fully illuminated (still 
neglecting the presence of spiders). The high-order loop is 
expected to run at 1 kHz. 

The AOF project undergoes in 2014-2015 a phase of lab¬ 
oratory tests with the aid of the dedicated ASSIST bench 
(Arsenault et al. 2013, 2014). The main idea is to test and 
to optimize the system before its delivery for the commis¬ 
sioning in Paranal. For this purpose, calibration and control 
tools have been developed in order to assist these tests with 
respect to e.g. DM/WFS mis-registration and turbulence 
diagnostics (Kolb et al. 2012; Garcia-Rissmann et al. 2013; 
Garcia-Rissmann et al. 2014), and are envisaged to be imple¬ 
mented as routines in the real-time computer (RTC) plat¬ 
form SPARTA (Suarez Valles et al. 2012). They are planned 
for use in either real-time or in post-processing modes, with 
the aim to improve the system performance by e.g. upload¬ 
ing optimized control matrices during operations (Kolb et al. 
2012 ). 

Through end-to-end system simulations, this work fo¬ 
cused on testing the reliability of a turbulence profiler in 
several aspects, which includes mimicking the observational 
conditions in Paranal. Section 2 describes the profiler code, 
how it operates and discusses a parameterisation used for 
the calibration of the profile. In section 3 we address some 
sensitivity tests for the assessment of the proHler perfor¬ 
mance, and section 4 describes its application to a more 
realistic scenario for all 3 modes. It is worth to point out 
that parameter dependencies such those on WFS flux and 
chosen centroid algorithms, although briefly studied as well 
(Garcia-Rissmann et al. 2014), are not addressed in this pa¬ 
per due to them being particularly inherent to the AO sys¬ 
tem performance, not to the profiler itself; this topic awaits 
a more dedicated study in the context of the overall qual¬ 
ity of the AOF correction. Section 5 summarizes and draws 
conclusions of the main results obtained in this work. 


2 THE Clr PROFILER 

Multi-WFS systems such as GRAAL and GALACSI were 
designed to probe the atmospheric turbulence in 4 different 
directions. The common volume of turbulence encompassed 
by pairs of these sensing directions makes the AOF - besides 
its main role in the AO correction - a suitable tool to trace 
the vertical turbulence distribution. The 40x40 SA^ array 
of each wavefront sensor, covering a pupil of 8 m across, 
can provide an unprecedented combined altitude resolution 
for this kind of application. In Open Loop (OL) regime, 
the 9920 measured slopes can be straightforwardly used to 
feed a turbulence prohler algorithm. In Closed Loop (CL) 
regime, Pseudo-Open Loop (POL) slopes for the profiler in¬ 
put can be obtained by multiplying the DM commands (c) 
by the interaction matrix (M) and adding this product to 
the residual slopes (sres): 

SPOL.fe = Sres.fc + M.Cfe_Afe (1) 



Figure 1. Left: AOF geometry for the four Shack-Hartmann 
WFSs. Solid and dotted lines indicate the high- and low-resolution 
baselines, respectively. The reference frame for x and y slopes di¬ 
rections is also shown. Right: a picture of one AOF WFS. Fully 
and partially illuminated SAs are denoted in white and grey, re¬ 
spectively. They altogether define the so-called “default” correla¬ 
tion mask, with 1240 SAs. 


where the subscript k refers to the frame number and Afc 
to the loop delay, i.e. the delay of the system between the 
application of the DM command and the availability of the 
new data from the WFS. This delay lies between 2 and 3 
frames for the AOF GLAO and LTAO modes. The interac¬ 
tion matrix M for these modes simply replicates the SCAO 
interaction matrix according to the number of WFSs. The 
following steps - valid for both OL and POL slopes - con¬ 
sist in a subtraction of the temporal averaged slopes and a 
frame-by-frame removal of the tip-tilt (TT) component from 
individual WFSs. 

Once the OL or POL slopes are correctly computed 
we can make use of some techniques to reconstruct the 
profile. These are post-processing tools that have been de¬ 
veloped and adapted from the conventional monitoring ones 
for the reconstruction of the turbulence distribution, based 
on AO RTC data. For instance, Wang, Schock & Chanan 
(2008) took advantage of the Palomar multiple guide star 
unit, cross-correlating its WFS data to obtain the distribu¬ 
tion of the Cn prohle as well as wind information. Later, 
Cortes et al. (2012) proposed two methods to obtain tur¬ 
bulent prohle distributions for the Gemini MCAO system 
GeMS (Rigaut et al. 2014; Neichel et al. 2014). The two al¬ 
gorithms are called ‘modihed-SLODAR’ and ‘wind-prohler’, 
the latter based on Wang et al’s work but adding the general¬ 
isation to laser beams. These methods have been re-baptized 
recently as ‘matrix inversion’ (MI) and ‘Fourier deconvolu¬ 
tion’ (FD) algorithms, respectively (Valenzuela et al. 2014). 

The aforementioned GeMS codes have been kindly pro¬ 
vided to ESO in the form of Matlab scripts; they have been 
adapted and modihed for the AOF case to account for its ge¬ 
ometry and to provide an alternative method for the unseen 
turbulence correction. After a phase of understanding the 
inner workings of the algorithms, only the wind-prohler one 
was kept, based on the robustness of the results obtained 
at the higher and the lower ends of the altitude range. This 
work is focused on testing this specihc algorithm, and we re¬ 
strict ourselves only to zero-delayed measurements (i.e. no 
wind speed estimation is attempted). For a recent compari¬ 
son of both methods applied to the AOF and a different ap¬ 
proach to the absolute calibration, see Valenzuela et al. 
(2014). 












Tests of a Cf Profiler for the AOF 3 


The AOF geometry is shown in Fig. 1, along with the 
directions we call x and y for the WFSs slopes reference 
frame. We can define two sets of resolution baselines: one 
refers to the larger angular separation connecting WFSO- 
WFSl and WFS2-WFS3 (shown with solid lines, that we 
call ‘high-resolution’ or HR baseline), and the other to the 
smaller angular separation connecting laterally the WFSs 
(shown with dotted lines, called ‘low-resolution’ or LR base¬ 
line). Fig. 1 (right) shows how the AOF WFSs look like 
in terms of valid (white-l-grey regions) and only fully illu¬ 
minated SAs (white region), defining correlation masks for 
these configurations (see explanation below). The grey pix¬ 
els denote the partially illuminated SAs, with fluxes above 
50 per cent of the those reached in fully illuminated SAs. 

The profiler algorithm relies on the cross-correlation be¬ 
tween WFSs that we call “A” and “B”. Its main equation 
is reproduced below (Guesalaga et al. 2014, and references 
therein): 


T^^{Au,Av,At) 





(t + At)'^ 


qAB 

'-^Au,Av 


( 2 ) 


where Su,v is the x or y slope contained in the SA denoted 
by the index (u, v) at time t, and Au and Av are displace¬ 
ments - in integer numbers of SAs - in the WFS grid. In this 
work, as mentioned previously, the frame delay denoted by 
At is always taken as zero. The summation ^ is carried 
out over all overlapping SAs, and the angle bracket symbols 
denote the time series average. The definition of the cor¬ 
relation masks is contained in Oau.Aui which compensates 
for the unequal number of correlated SAs as one moves away 
from the center of the correlation maps. For the AOF, which 
suffers no fratricide effect, and in particular for this set of 
simulations without the presence of spiders, is the same 
for all pairs of WFSs^ 

The next step comprises a deconvolution for the impulse 
response of one single average WFS: 

TAB* = jr-l ^ Q gjj ^3^ 

where T and T~^ denote the Fourier and the inverse Fourier 
transforms, respectively, and jg deconvolved cross¬ 
correlation map. The profile is extracted from the decon¬ 
volved maps along the directions that connect the baselines, 
starting from the center. The redundant baselines sharing 
the same separating angles of laser guide stars are combined 
to provide the HR and LR measured profiles. 

The absolute calibration of the above HR and LR pro¬ 
files requires reference data, which is accomplished by sim¬ 
ulating a large number of independent phase screens with 
pre-defined turbulent parameters. The WFS slopes derived 
from such phase screens are combined to emulate a flat pro¬ 
file observed in a large number of frames, from which the 
average angle-of-arrival (AA) variance is derived. A similar 
analysis as that carried out for the measured data is per¬ 
formed, obtaining a calibration profile which can be used 
to normalize the measured one. Finally, multiplying such 


^ This assumption has to be re-evaluated for the ASSIST bench 
or telescope measurements because of the presence of spiders. 


normalized measured profile by the AA-variance from the 
calibration data provides the stratified distribution of slope 
variances in the atmosphere. This is then corrected by the 
cone effect and converted into a distribution of Fried pa¬ 
rameters per altitude bin using a parameterisation of the 
von Karman atmospheric model of the type: 

criA oc X /(Lo) (4) 

in units of rad^, with d = 0.2 m (SA size), A = 589 nm and / 
being a function of Lo only. Notice that the calibration pro¬ 
cedure makes an assumption on the value of Lo that should 
reflect the outer scale of the analysed data. The distribution 
of outer scales measured in Paranal shows a concentration 
around 22 m (Martin et al. 2000). The AA-variance from 
a von Karman model has a strong dependence on this pa¬ 
rameter on smaller scales, compromising the absolute profile 
calibration if a good estimate of Lq from data is precluded. 
For this work we used a parameterisation derived by Conan 
et al. (2000) - with / presenting a dependence on Lq given 
as powers of -1/3, -2 and -7/3, and checked against the AA- 
variances measured from the calibration phase screens - in 
contrast to the original work on GeMS, which adopted a 
Kolmogorov model (/=1). 

The ro distribution obtained from the adopted param¬ 
eterisation is easily converted to the HR and LR profiles 
of the seen turbulence, in units of rrd^^. On the other hand, 
the integration of these profiles and their conversion to 
A = 500 nm can provide the relevant Fried parameter of the 
integrated seen turbulence for both resolutions. In order to 
get rid of any edge effects that occur in the cross-correlation 
maps, only bins from 1 to Nhr and Nlr are kept, dis¬ 
carding those of higher altitudes. For the GLAO modes and 
assuming that the analysis was carried out on all valid SAs 
(i.e. using the default mask), as default we build a combined 
profile by keeping the first 38 HR bins and populating the 
higher altitudes with the remaining selected LR bins. This 
generates 46 altitude bins in total, from which a combined 
Fried parameter can also be derived. For the LTAO mode we 
discard the LR profile - too poor in resolution - and keep 
only the 20 first bins of the HR one. The main parameters 
for the profiler in each mode are shown in Table 1. 

The altitude bins (in metres) can be calculated from the 
following formula (for both HR and LR baselines): 


hm 


m d h-Na 

sec 2 h-Na 6' + m d 


+ ho 


m = 0, l...Nb — 1 


( 5 ) 


where sec 2 corresponds to the airmass of the observation 
(the secant of its zenith angle z), /iNa to the altitude of the 
Sodium layer, 8' to an angle - in radians - related to the 
baseline under consideration (from Table 1, we set 0' = Shr 
and 9' = 8\,^l\f2 for HR and LR baselines, respectively). 
The parameter d (=0.20 m) is again the SA size and Nb 
the number of sampled bins (40 for HR and 28 for LR). 
Finally, the parameter ho is a zero point scale for the binned 
altitudes. In the case of the AOF, ho = —89 m accounts for 
the conjugation of the entrance pupil to this altitude below 
the ground, which is associated to the DSM location. 

Fig. 2 shows hmax, defined as the maximum height 
achieved by the profile, as a function of the zenith angle, 
for only the AOF GLAO modes (GRAAL and GALAGSI 
WFM). With the two GALAGSI modes we are able to cover 
- although less strictly in WFM - the bulk of the turbulence 
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Table 1. Maximum altitudes and spatial resolutions probed by the profiler using all 1240 valid SAs (default mask). 
NFM: Narrow Field Mode; WFM: Wide Field Mode. Columns 3 & 7: angular separation of baselines. Columns 4 & 
8: number of adopted bins for baselines. Columns 5 &; 9: maximum altitudes for baselines given by the profiler, at 
2 = 30°. Columns 6 & 10: ranges of spatial resolution in baselines, at 2 = 30°. 





Low 

Resolution 




High Resolution 



Mode 

^LR 

Nlr 

lz=30° 

^max,LR 


6hr 

Ahr 

i,z=30° 

^max,HR 






[km] 

[km] 



[km] 

[km] 

GALACSI NFM 

LTAO 

14.1” 

27/0“ 

45.5/-“ 

3.4-0.9/-“ 

20.0” 

38/20' 

* 37.9/24.5“ 

1.7-0.6/1.7-0.9“ 

GALACSI WFM 

GLAO 

1.51’ 

27 

12.4 

0.55-0.41 

2.13’ 

38 

9.14 

0.28-0.22 

GRAAL 

GLAO 

8.20’ 

27 

2.49 

0.102-0.096 

11.6’ 

38 

1.76 

0.051-0.049 


“The slash symbol separates the maximum possible and the actually adopted number of bins, altitudes or resolution 
ranges for GALACSI NFM, respectively. In this case, only the HR bins are considered in the analysis. 


AOF GLAO SYSTEMS 



Figure 2. Maximum altitudes probed by the AOF GLAO modes 
as a function of the zenith angle, assuming that all valid SAs 
(black) and only the fully illuminated ones (grey) are used for 
the correlation mask. The AOF LTAO mode is not shown here 
since its maximum altitude reaches much beyond 20 km for both 
mask cases, no matter the zenith angle of the observation (0° Si 
2 60°). 


contained in the troposphere. GRAAL offers extremely good 
spatial resolutions (as shown by (5hr and 5lr for 2 = 30°, 
see Table 1) but fails to provide information for altitudes 
higher than 2.9 km (or ~1.4 km, in the extreme case of 
2 = 60°). As mentioned previously, with GALACSI NFM 
we limit ourselves to 20 HR bins for them being enough to 
probe, with the best resolution, up to 24.5 km at 2 = 30°. 
This way we avoid degrading our results with the noisier 
and negligible higher turbulence layers. 

The simulations of all ‘measurement’ data shown in this 
work have been generated with the Octopus end-to-end sim¬ 
ulation tool. This is a parallelized code that runs on a Linux 
cluster dedicated to simulations of adaptive optics systems 
at ESO. A good comparison of Octopus with other existing 
simulation tools in the AO community can be found in Le 
Louarn (2010). Throughout this work we adopted a pupil 
size of 400 pixels to represent the 8 m telescope aperture. 


with a central obscuration of 14.5 per cent. The AOF WFSs 
contain SAs with 6x6 pixels, with a field of view (FoV) of 
0.81” pixel”^. The Low Light Level CCD detector was simu¬ 
lated by doubling the amount of photon noise and adding 0.2 
electrons per pixel of background noise to take into account 
the readout noise left after the amplification of the detector 
signal. For simulations mimicking realistic conditions (sec¬ 
tion 4) 10 turbulent layers were adopted, distributed up to 
14 km and allowed to displace in random orthogonal direc¬ 
tions with user-defined turbulence weights and wind speeds. 
Turbulent phases have been generated on a support of 8192 
pixels length, meaning that they wrap up for wind speeds 
larger than 11 ms”'^ in typical simulation timespans (15 
seconds). However, because such layers are concentrated in 
higher altitudes and superpose only with smaller weights to 
the total turbulence, any repetition effects should be negligi¬ 
ble. For test cases such as the single layer analysis (see next 
section) the wrapping of phase screens would, at maximum, 
set a conservative error to the results. 


3 PERFORMANCE ASSESSMENT CASES 

GALACSI WFM, the intermediate spatial resolution case, 
is the mode under inspection in this section. Other AOF 
modes are treated later, once we have an insight on the lim¬ 
itations of the profiler in this intermediate case. In order to 
disentangle the effects introduced by the POL computation 
from GLAO/LTAO CL modes and the outcome from the 
profiler, we opted for analyzing in this section only OL sim¬ 
ulations. Taking this approach, we bear in mind that we may 
incur some non-linear behaviour of the WFSs. All simula¬ 
tions in this section have been run in the WFS high photon 
flux regime. Turbulence parameters of (ro,Lo) = (0.10, 25) 
m have been used. We also assumed 2 = 0° so the probed 
altitudes were at their maximum, and /lo = 0 m. 

3.1 Response to Different Turbulence 
Distributions 

3.1.1 Single-Layer Profiles 

For this section, simulations comprised of profiles of one sin¬ 
gle atmospheric layer concentrating all turbulence have been 
analysed. The layers were simulated at the exact altitudes of 
the GALACSI WFM HR bins as given by equation (5) with 
O' — 0hr. Notice that turbulent layers beyond the altitude 
corresponding to Nb = 56 are not supposed to be probed 
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HR Layer #9 @ 2506.5 m 


HR Layer #48 @ 12965 m 




HR Layer #49 @ 13201 m 



HR Layer #65 @ 16781 m 



Figure 3. Examples of the results from the profiler for simulations of single layers at different altitudes. The grey bar represents the 
input value (=100 per cent in one layer) whereas the crosses and dots show the output results as given by the profiler and noise-filtered, 
respectively. Layer #9: exact HR bin; layer #48: falling between two high altitude LR sampled bins; layer #49: falling at an exact 
sampled high altitude LR bin; layer #65: altitude not sampled by the profiler. The apparent high levels of noise, produced by a direct 
normalization of the profile, can be corrected by taking into account the unseen turbulence (see later). 


by the profiler. Fig. 3 shows typical examples of the profiler 
response when having the turbulent layer located at the 9**', 
48*^*', 49*’’ and 65*’’ HR altitude bins, with the integrated 
output from the profiler normalized to 1. Input layers 
are given by the grey bars and the raw output from the 
profiler by the crosses. In terms of qualitative response, the 
profiler succeeds in reconstructing the input layers at their 
locations. We notice different types of response: 


(i) Strong response in one output layer when the input 
layer falls exactly in a HR or a LR profile bin (e.g. layers 
#9 and #49, respectively). 

(ii) Strong response in 2 layers when the input layer falls 
in a HR-like bin not sampled by the HR part of the pro¬ 
filer but it is still within the range of the combined profile 
(e.g. layer #48). In other words, the sum of adjacent bins 
approximately matches that of the input profile. 

(hi) No response when the input layer is beyond the re¬ 
constructed profile (e.g. layer #65). In this case the normal¬ 
ization of the raw output profile to 1 explains the apparent 
high level of noise seen in the figure. The noise, in abso¬ 
lute Cn value, is actually the same as in the other types 
of response. Once the unseen turbulence is accounted for 
(see discussion later), the apparent noise in the normalized 
profile drops significantly, looking the same as in the other 
cases. 


Noise is observed in all layers where no turbulence is ex¬ 
pected, resulting in an output value at the location of 
the peak with less than 100 per cent. By selecting a certain 
number of bin data points that would exclude the turbulent 
ones (for instance, the N ^ Nt — IQ smallest values of the 
output profile, where Nt is the total number of combined 
bins probed by the profiler), we can compute their median 
Mp = median(P) and median deviation values defined by 

crp = \/ median(P — MpY (6) 

and devise a filter that sets to zero all values lower than 
Mp -f a.ap. The parameter o is a threshold to be defined 
and that depends on the amount of noise to be filtered. The 
purpose of this procedure is to remove unrealistic profile 
values e.g. those with negative intensities. Once the filtering 
is applied, the should be renormalized to 1. Applying the 
filter (with a = 1) to every profile obtained from simulated 
data with one layer at a visible altitude (<14 km), we obtain 
median values of the noise which are very small (0.13 per 
cent on average) and oscillating close to 0 per cent. The 
median deviation of the noise varies between 0.4 and 0.6 per 
cent with an average value of about 0.5 per cent. This value 
can be viewed as the noise ‘standard deviation’ of the 
profile estimation in layers without turbulence. On the other 
hand, the ratio of the output turbulence in a determined 
layer to the total turbulence - estimated from a simulation 
with 100 per cent input turbulence exactly at this bin (or 
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in the 2 adjacent layers when in LR) - lies between 75- 
90 per cent for the HR bins and around 100 per cent for 
the LR bins. In other words, the peaks are detected in all 
resolutions but their strength is lower than the input one, 
specially in HR. A correction with a a = 1 Hlter can make 
the results worse. This shows that by not hltering enough 
the noise the whole profile is biased towards positive values, 
decreasing the relative amplitude of the peak; for a noise 
Hlter with a = 2 the results turned out to be similar to the 
non-filtered ones. 

Bear in mind that the layers above ~14 km are not sam¬ 
pled. One should filter even more the noise or broaden the 
bins but that may not be a valid solution when using real 
data for losing resolution. The better response of the higher 
layers suggests that the LR profile is less noisy and gives 
more accurate response to the input turbulence. A sugges¬ 
tion would be that, when accuracy is more important 
than spatial resolution, one should use the LR profile only; 
this could be the case, for instance, of GRAAL profiles. 

If the turbulence is distributed above and below the 
highest probed altitude hmax, the Fried parameter de¬ 
fined as the one measured by the profiler, is overestimated. 
On the other hand, estimates made from the WFSs AA- 
variances - properly corrected by the cone effect - should 
reflect the total turbulence along their lines-of-sight, being 
closer (ideally equal) to the input value of the simulations^. 
Equation (4) for a von Karman model implies that the com¬ 
putation of these parameters depend on an assumption for 
the I/O value of the input data, a tricky parameter to both 
simulate and estimate with the current algorithm. This has 
an impact on the absolute calibration of the profile but for¬ 
tunately not on its shape. For that reason, and because most 
applications of the profiler to optimize the AOF performance 
only require information about the relative distribution, 
m this work we restrict ourselves only to the analysis of the 
normalized turbulence profiles. This means that and 
in this work are just used for a correction of the unseen tur¬ 
bulence. We will return to this discussion when analyzing 
multi-layer profiles. 


3.1.2 Multilayer Profiles 

A natural extension to the single layer analysis is the one 
of a few adjacent layers. This is useful to gain an insight on 
the cross-talk between layers in the profiler response. There 
is some evidence of this behaviour in Fig. 3, with turbu¬ 
lence input in bin where we see a small ‘leakage’ of the 
layer strength to the neighbouring higher altitude bin. In¬ 
put profiles composed of 1, 2, 3, 4 or 5 layers (containing 
each an equal share of the total turbulence) located at the 
exact altitudes of an output HR profile bin close to 5 km 
and the bins immediately above have been simulated to fur¬ 
ther check this behaviour. With a = 0 for the noise filtering, 
the integrated value of the output turbulence in those layers 
tends to increase from 73 to a maximum of ~84-85 per cent 
with the increase in the number of turbulent bins from 1 to 
5 (Fig. 4). The results also show some leakage of intensity to 

^ This still might not be true because of the limited time sampling 
and linearity/spot truncation effects of the centroid algorithms. 
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Figure 4. Top: integral of the normalized profile, only in 
the turbulent adjacent bins. Bottom: integral of the normalized 
profile, in the vicinity of the turbulent bins. Curves show the 
results for an applied noise filter of « = 0, 1 and 2. The results 
are for the HR part of the profile. 


neighbouring layers that can reach ~7 per cent, tending to 
decrease for a larger number of input turbulent bins. Using 
a = 1 and 2 for the noise filtering produces slightly better 
results, although with final integrated intensities still below 
90 per cent. In summary, this test has shown that the pro¬ 
filer seems to respond better to a turbulence which is not 
concentrated in isolated thin slabs, at least in the HR part of 
the profile (for LR the convergence is better, as mentioned 
earlier in the single-layer case). 

A case of a fiat input profile composed of Nb = 80 layers 
has also been analysed, with each layer located at the ex¬ 
act altitudes of the HR profile bins. In practical terms, this 
means that there is a significant turbulence fraction which 
is not seen by the profiler. A first inspection shows that the 
average level of the estimation turns out to be good, but 
the output profile seems to increase between ~10-14 km. 
This happens because at those altitudes each output bin 
(which is LR) covers 2 layers of the input profile thus dou¬ 
bling their total turbulent strength. A much fairer represen¬ 
tation consists in resampling the simulation layers intensi¬ 
ties to the output bins, assuming a linear interpolation based 
on results from tests with sub-bin displacements. To do so, 
each simulated layer intensity is spread into the 2 nearest 
output bins i and i + 1 with a ratio proportional to the 
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Figure 5. Response of the profiler to a flat distribution. Top: 
uncorrected by unseen turbulence; bottom: corrected by unseen 
turbulence, using the prior knowledge of the distribution from 
the simulations. The rms between the output profiles minus the 
resampled input ones is 1.12 and 0.7 per cent for top and bottom 
cases, respectively. 

distance of hi+i and hi to that layer altitude, respectively. 
For instance, if Hk (= input layer altitude) lies between hi 
and /i 2 (i.e. hi ^ Hk ^ h 2 ), the intensity of the input bin 
will be shared between the two output bins with weights of 
{h 2 —Hk)/{h 2 — hi) and {Hk — hi)/{h 2 — hi), respectively. In 
this representation, input contributions which lie above the 
maximum probed altitude remain unaffected. Fig. 5 (top) 
shows the resampled input and the output profiles. The rms 
of their difference in the sampled bins is of 1.12 per cent 
(almost the intensity of the input signal in each of the 80 
HR bins, i.e. of 100 / 80 ~ 1.30 per cent). 

A still better match of the output profile can be ob¬ 
tained if we compensate the output profile by the unseen 
turbulence. In this specific case, the Fried parameter input 
in the Octopus simulations, called hereafter ro™, is 10.3 cm. 
On the other hand, is 14.6 cm and the average Fried 
parameter obtained from the WFSs AA-variances, still un¬ 
corrected by the cone effect, is 12.5 cm. An equivalent tur¬ 
bulence altitude can be estimated using the formula: 

( 7 ) 

where the input distribution can come either from the 
profiler output or from a priori knowledge of the input 
distribution (e.g. from simulations). The cone effect factor 



Figure 6. Hypothetical profile used in the analysis from sections 
3.2 to 3.4, resampled to the output bins of the profiler. The output 
profile shown here refers to the profiler results for a GALACSI 
WFM OL simulation of high flux (800 photons SA“l frame“'^), 
an undersampling of t =10 and a dataset length of 15000 frames. 
A noise filtering of o = 1 was applied. 

can then be approximated by /c ~ 1 — sec« h/h^a,. Such a 
factor produces a estimate of 11.3 and 10.9 cm, if we 
consider h equal to 8.47 km (using all 46 output profile bins 
and intensities) and 11.8 km (using all 80 seen-|-unseen simu¬ 
lated bins and intensities) for such a correction, respectively. 
In the case of real measurements, the output profile should 
thus be multiplied by to account for the un¬ 

seen turbulence. Notice at this point that any dependence on 
Lq is cancelled out in a normalized profile. Fig. 5 (bottom) 
shows the improvement obtained in the output profile if we 
apply the correction considering all simulated bins. The rms 
of the residuals after the cone effect and unseen turbulence 
correction drops to 0.7 per cent in this case. Naturally the 
improvement can be even better if one provides an indepen¬ 
dent and better estimate of the total ro, since the normal¬ 
ization factor of (14.6/10.3)“®^® = 0.75 leads to a better 
correction of the unseen turbulence (rms of the residuals ~ 
0.5 per cent). In summary, even not constituting a realistic 
scenario because of its profile flatness, this case is helpful 
to validate the bin resampling procedure and the correction 
for the unseen turbulence. It confirms very clearly the 0.5 
per cent figure for the noise standard deviation of the profile 
estimation in this extreme case of large equivalent altitude 
turbulence. It is also worth pointing out that the only AOF 
mode likely to be significantly affected by unseen turbulence 
is GRAAL because of its low altitude coverage (maximum 
height of 2.9 km at zenith). 

3.2 Dependence on Temporal Data Sampling 

The quantities used in the derivation of the profiles are 
temporal statistical averages. Therefore, it may be possible 
to achieve in the post-processing stage the same results using 
only a smaller subset of the original number of frames. 

A hypothetical profile as shown in Fig. 6 is adopted for 
this test. As a quality criterium we use the median abso¬ 
lute deviations of the residuals (resampled normalized input 
profile minus the normalized output profile). Since this hy- 
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the output when the number of acquired frames is less than 
~10000 (@ 1 kHz loop rate). For comparison, the GeMS 
circular buffer files analysed in Cortes et al. (2012) had a 
maximum size capacity of 24k frames which, in a loop rate 
tuned to compensate for the photon return (<800 Hz), pro¬ 
duced datasets of the order of a few minutes. Regarding the 
CANARY bench working in MOAO, Morris et al. (2013) re¬ 
ported an acquisition of 2-3 minutes for the determination of 
the profile from WFSs telemetry data. Being conservative, 
we adopt 15 seconds as our minimum individual ‘acquisition’ 
time at the nominal loop rate. In real observations, a pro¬ 
cedure of averaging profiles obtained from different datasets 
spanning a couple of minutes might be adopted as well in 
order to ensure a better statistical result. 
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Figure 7. Statistics of the residuals: median and maximum abso¬ 
lute deviations for different samplings. Top: taking only a certain 
percentage of frames in the cross-correlation computation; bot¬ 
tom: taking different dataset lengths, for a fixed percentage of 
frames (r = 10). 


3.3 Dependence on Input Lo 

The Octopus simulation studied in this case is the same as 
the one in the previous section. The analysis, however, dif¬ 
fers in the assumption of the Lo for the computation of the 
profile, ranging from 5 to 100 m. The normalized output 
profile is not sensitive to the assumed calibration outer scale, 
as one could expect, being consistent with median absolute 
deviations in the residuals of 0.5 per cent. The sensitivity to 
this parameter is pronounced in the estimation of and 
reaching deviations of up to ~35 per cent when the as¬ 
sumed I/O is very different from the effective Lq of the simu¬ 
lations. This confirms that this profiler method alone should 
not be used to estimate the Fried parameter. An absolute 
calibration of the output profile through this method would 
then require an independent estimate of turbulence param¬ 
eters. Notice, however, that a newer and smart approach 
to calibrate the absolute profile has been addressed in 
Valenzuela et al. (2014); it is based on fitting the covariances 
profiles, being independent of prior model assumptions. 


pothetical profile has no turbulence contribution beyond 14 
km, no unseen turbulence correction is needed, and there¬ 
fore the integrated turbulence in the probed range is 100 
per cent. For the case shown in the figure, considering only 
bins with non-zero input turbulence, the adopted residual 
criterium gives 0.9 per cent. 

Fig. 7 (top) shows our statistical criterium as a function 
of the adopted undersampling interval for a dataset length 
of 15000 frames (15 seconds @ 1 kHz). An undersampling 
of r means that only 1/r of the total number of frames are 
taken into account in the cross-correlation computation. It 
can be clearly seen that for up to an undersampling interval 
of about 10-50 frames the residuals do not change much but 
start to do significantly for intervals higher than about 100 
frames. The undersampling parameter has an impact on the 
time required to compute the profile and this can be of great 
value for the real time optimization of the AO correction. 

Another possible test refers to the dataset length. This 
aspect has an impact on the minimum time required for the 
acquisition of telemetry data in order to retrieve a reliable 
atmospheric profile estimate. Fig. 7 (bottom) shows the re¬ 
sults for this test, with dataset lengths ranging from 6000 to 
30000 frames and adopting a fixed undersampling interval 
of 10 frames. Apparently the dataset length starts to affect 


3.4 Dependence on Input Na Layer Altitude 

On the real system, the RTC provides information on the 
average Sodium layer altitude through the trombone and fo¬ 
cus compensator loop. As an input to the profiler, the same 
information impacts on the definition of the altitude bins 
through equation (5). It has been verified - with the same 
simulation data of Fig. 6 - that when an error is made on 
the assumed altitude of the Na layer, the positions of the 
layers are mistaken by a small fraction of that amount and 
the integrated profile amplitude changes slightly. For exam¬ 
ple, an error of 10 per cent on the input Na layer altitude 
leads to an error of 1.4 per cent on the layers position and 
to a maximum absolute deviation of 5 per cent from the ex¬ 
pected amplitudes. This also implies that the estimation of 
and are little affected as well, because their depen¬ 
dencies on the integrated profile and on the cone effect cor¬ 
rection variations are negligible. In summary, the profiler 
does not require a good accuracy on the Na layer distance; 
for this estimation one could assume a median altitude of 
90 km and the actual zenith angle, with no need of informa¬ 
tion from the actual position of the trombones or the focus 
compensators in the AO modules. Notice, however, that this 
information does matter for the optimal system correction. 
Any uncertainties in the Na layer altitude during operations 
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POL Slopes Variance - GALACSI NFM 



Figure 8. Left: Example of POL slopes variances - in pixel^ 
units - measured in a WFS from simulations of the GALACSI 
NFM, under a typical Paranal turbulence of 0.6” seeing. Right: 
in white, the mask designed to cut off problematic correction 
regions, containing 1044 SAs to be used in the profiler. 


can indirectly affect the estimation because of the loss of 
focus in the WFS cameras and its consequent implications 
for the slopes estimation (introduction of noise, truncation 
effects, etc). 


4 PERFORMANCE UNDER TYPICAL 

PARANAL CONDITIONS 

Most of the tests shown in previous sections have been car¬ 
ried out using GALACSI WFM in OL, simulated at zero 
zenith angle and taking into account the conjugation of the 
telescope pupil to the ground. Hereafter more realistic as¬ 
sumptions are adopted, evaluating CL performances under 
typical telescope setup and observatory turbulence condi¬ 
tions. The actual value for the altitude conjugation of the 
DSM is then set to ho — -89 m. The flux on the WFSs 
is taken as the nominal one for the AOF, of 80 photons 
SA“^ frame”^. This is a safe photon flux regime for robust 
profiler results as shown by previous tests (Garcia-Rissmann 
et al. 2014). The zonal interaction matrix used in equation 
(1) for computing the POL slopes is calculated through Oc¬ 
topus in a SCAO configuration and replicated to account 
for all 4 WFSs when transforming commands into slopes. In 
GALACSI NFM, ground and ‘virtual’ commands provided 
by the Octopus LTAO simulations are summed to give one 
single set of commands to be multiplied by M. It is impor¬ 
tant to point out that the reconstructor and the loop gain 
have yet to be fine-tuned for the maximum Strehl, but sta¬ 
bility in the loop was ensured. 

Fig. 8 shows the variance of POL slopes mapped to 
one of the WFSs, as seen in a CL LTAO case for a typi¬ 
cal Paranal profile with seeing 0.6” (@ 0.5 pm). The figure 
shows an increased variability in the WFS edges, reflect¬ 
ing a larger uncertainty of the correction in those regions, 
partially because of discontinued actuators dynamics at the 
borders. Notice that a wavefront aberration improperly cor¬ 
rected by the DSM at the pupil edges leaves a footprint com¬ 
mon to all sensing directions, and this feature translates into 
a strong cross-correlation when there is an overlap of WFSs 
maps (i.e. for the lowest altitude bins). Some preliminary 
tests have shown that even the fully illuminated SAs mask 
is not effective in dealing with such high variance regions. 


Table 2. List of typical turbulence profiles compiled in Paranal. 
/3: scaling factor for the layer speeds; Prob: probability of occur¬ 
ring such conditions in Paranal; T.Q.: turbulence quality. Notes: 
ro @ 0.5 pm; tq @ 2 : = 30'’. 


ID 

S 

ro 

TO 

P 

h 

Prob. 

T.Q. 


[”] 

[m] 

[ms] 


[km] 

[%] 


POl 

0.4 

0.186 

4.6 

1.29 

3.73 

7.0 

median 

P02 

0.6 

0.136 

3.9 

1.44 

2.39 

6.0 

good 

P03 

0.6 

0.136 

3.8 

1.14 

3.73 

12.0 

median 

P04 

0.6 

0.136 

3.9 

0.92 

4.88 

6.0 

bad 

P05 

0.8 

0.116 

3.0 

1.44 

2.66 

6.5 

good 

P06 

0.8 

0.116 

3.0 

1.15 

3.82 

13.0 

median 

P07 

0.8 

0.116 

3.1 

1.00 

4.45 

6.5 

bad 

P08 

1.0 

0.101 

2.4 

1.44 

3.05 

4.5 

good 

P09 

1.0 

0.101 

2.5 

1.19 

3.81 

9.0 

median 

PIO 

1.0 

0.101 

2.4 

1.12 

4.47 

4.5 

bad 

Pll 

1.2 

0.089 

2.0 

1.55 

2.73 

3.0 

good 

P12 

1.2 

0.089 

2.1 

1.27 

3.53 

6.0 

median 

P13 

1.2 

0.089 

2.0 

1.36 

3.69 

3.0 

bad 

P14 

1.4 

0.074 

1.4 

1.64 

3.53 

13.0 

median 


Therefore, in order to circumvent these spurious effects a 
new mask based on the variance of the slopes - see also 
Cortes et al. (2012) - is adopted, with 1044 ‘good’ SAs (Fig. 
8, right). In practical terms, masking out the problematic 
SAs restricts the altitudes probed by the profiler, shorten¬ 
ing them by about 1 bin in maximum height (LR high end 
of the profile). Applying this mask to GALACSI NFM POL 
slopes has no limiting effect on its output HR profiles be¬ 
cause of the large altitude range achieved in this mode. For 
GLAO modes the maximum height loss is of about 100 m 
and 400 m for GRAAL and GALACSI WFM, respectively. 

Once defined the proper setup, we proceed to present 
the realistic profiles to be fed to the profiler in the CL tests. 
Table 2 shows a compilation of typical profiles obtained 
in Paranal from a mix of SLODAR, MASS-DIMM and SCI- 
DAR measurements collected along several years (Sarazin 
et al. 2013). Groups have been divided according to seeing 
bins in terms of ‘median’, ‘bad’ and ‘good’ profiles, based 
on their turbulence contributions closer to the ground, i.e. 
those easier to be corrected by the AOF. These groups have, 
respectively, median values of h (equation 7) of 3.73, 4.46 
and 2.70 km. The 7*’*' column lists the probability of observ¬ 
ing such profiles in this database. In the simulations both 
the turbulent heights and the Na layer altitudes/thicknesses 
(with default values of 90/7 km, respectively) are multiplied 
by sec 2 , given that we assume here a default zenith angle 
of 30°. Also as input for the simulations, the reference wind 
profile Vref has to be scaled by a constant /3, given in Tables 
2 and 3, such that 


f EJhVref(h)]^/^C^(h)Y^^ 

[ E,c^(h) J 


( 8 ) 


where the relation 

TO = 0.314 ^ (9) 

V 

holds. For the seeing S and coherence times ro shown in 
Table 2 the P values have been calculated and used in the 
Octopus simulations. 

Finally, we input an outer scale of 25 m in all simu¬ 
lations and restrict ourselves to analyzing only normalized 
profiles for the reasons exposed previously. 
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Table 3. Typical Paranal turbulence profiles used in the simulations. Pn: profile identification, as in table 2. 


Altitude 

'^ref 

POl 

P02 

P03 

P04 

P05 

P06 

P07 

P08 

P09 

PIO 

Pll 

P12 

P13 

P14 

m 

[ms-i] 

[%] 

[%] 

[%] 

[%] 

[%] 

[%] 

[%] 

[%] 

[%] 

[%] 

[%] 

[%] 

[%] 

[%] 

30 

5.7 

70 

83 

70 

53 

77 

59 

41 

65 

45 

16 

46 

26 

0 

26 

140 

5.1 

1 

1 

1 

1 

2 

2 

1 

5 

4 

2 

10 

8 

3 

8 

281 

4.4 

3 

3 

3 

2 

4 

4 

4 

7 

7 

5 

12 

11 

6 

11 

562 

3.9 

5 

5 

5 

5 

5 

6 

8 

7 

9 

12 

15 

16 

23 

16 

1125 

4.4 

0 

0 

0 

0 

0 

1 

1 

0 

3 

11 

1 

6 

28 

6 

2250 

7.2 

2 

0 

2 

7 

1 

5 

13 

1 

6 

21 

1 

10 

14 

10 

4500 

14.2 

5 

1 

5 

11 

3 

9 

15 

4 

12 

16 

6 

10 

14 

10 

7750 

26.3 

4 

2 

4 

6 

2 

4 

5 

4 

5 

6 

4 

6 

5 

6 

11000 

29.8 

4 

2 

4 

6 

3 

5 

6 

4 

5 

6 

3 

4 

4 

4 

14000 

15.2 

6 

3 

6 

9 

3 

5 

6 

3 

4 

5 

2 

3 

3 

3 
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Figure 9. Upper plots: performance of the profiler for GRAAL with the 14 typical Paranal profiles in OL (left) and CL (right). Solid 
lines: median absolute deviations; dashed lines: maximum absolute deviations; dots: o = 0 noise filtering; stars: a = 1 noise filtering; 
open circles: a = 2 noise filtering. Thick lines in the top-right plot refer to results obtained using the mask of only 1044 SAs and 0 = 1 
noise filtering. Lower plots: two examples of outputs from the profiler using CL simulations: of a good profile with 0.6” seeing and 8 per 
cent of unseen turbulence (left) and of a bad profile with 1.2” seeing and 26 per cent of unseen turbulence (right). All output profiles have 
been corrected by the unseen turbulence using as a reference the input value of the simulations, rg™. Noise filtering of o = 1 applied to 
the output profiles in the lower plots. 


4.1 GLAO Modes 


4.1.1 GRAAL 

GRAAL is one of the GLAO or ‘seeing-enhancer’ modes 
for the AOF. It is aimed to feed the 7.5’x 7.5’ FoV imager 
HAWK-I on the Nasmyth A focus of UT4. Its goals are an 


improvement of about 40 per cent in the A-band FWHM and 
to achieve observations with 0.3” F’WHM in ~50 per cent 
of the time under a seeing slightly below 1” (Paufique et al. 
2012). As shown in Table 1, the separation angle for the HR 
and LR baselines are 11.6’ and 8.2’, respectively. This gives 
the profiler an unprecedented resolution in probed altitudes 
but a maximum altitude constrained to less than 2.89 km 
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at zenith, and 2.49 km at z = 30°. Therefore, for the typical 
Paranal listed profiles with equivalent altitudes h ranging 
between 2.39 and 4.88 km, the amount of unseen turbulence 
can be significant, lying between 8-33 per cent for 2 = 30°. 

Fig. 9 (top) shows the performance of the profiler for all 
the normalized profiles listed in Tables 2 and 3, under differ¬ 
ent noise filters (a = 0, 1 and 2). For the sake of comparison, 
the results for profiles derived from OL simulations are also 
shown on the left. The statistics are carried out only on 
resampled input bins that contain turbulence. Median and 
maximum absolute deviations lie between ~0.5-2.5 and ~5- 
20 per cent, respectively, and do not depend significantly 
on the noise filter applied. There is a general trend for the 
error in the peak estimation to decrease for larger tur¬ 
bulences, in similar ways for profiles derived from either OL 
or CL simulations. Subsets of profiles with the same seeing 
tend to show a better performance (smaller maximum resid¬ 
uals) for a smaller turbulence content close to the ground, at 
least for seeing conditions better than 1.2”. This is related 
to the trend to underestimate turbulence peak intensities, as 
discussed earlier. Running the profiler with the mask shown 
in Fig. 8 produces no significant effect on the statistics of 
the residuals, as shown by the thick lines in the CL case. 

Fig. 9 (bottom) shows two opposite examples of profiles 
retrieved from CL GRAAL data, a good turbulence profile 
with 0.6” seeing and 8 per cent of unseen turbulence, and 
a bad turbulence profile with 1.2” seeing and 26 per cent 
of unseen turbulence. The normalized output profiles were 
noise filtered with a = 1 and unseen turbulence corrected us¬ 
ing Tq™ as a reference; they have been retrieved using both 
the default and the 1044 SAs masks (denoted by crosses 
and open circles, respectively). The very two first bins of 
the GRAAL profile refer to altitudes conjugated below the 
ground (-89 and -38 m) but we see that, using the default 
mask, the contribution obtained in those bins are not neg¬ 
ligible, amounting to ~ 12-23 per cent in the shown cases. 
It is important to highlight that this effect is not seen in 
OL data, and therefore it is related to the DSM correction 
and/or to the POL reconstruction. By applying the mask of 
1044 SAs we managed to reduce such amount to ~6-16 per 
cent, with the contribution in the second altitude bin staying 
practically unaltered. This might indicate that either a more 
restrictive mask is needed and/or that there is a cross-talk 
effect between turbulent bins, as seen previously. In general, 
despite the detected biases in the negative altitude bins we 
see that the profiler succeeds qualitatively in identifying the 
major turbulence peaks and the overall shape of the input 
profile. The underestimation of the output CL profile (using 
the default mask) in specific peaks can be improved if we 
set as zero the turbulence in the two first bins and re-do the 
normalization. For instance, in the 0.6” case, the peaks at 
altitudes 13 m and 62 m reduce their residuals from ~10 
and 20 to -3 and 18 per cent, respectively. Given all the un¬ 
certainties, errors in the retrieved unseen turbulences for all 
the profiles can reach 16 per cent, with an average value of 
8 per cent. 

4.1.2 GALACSI WFM 

The other AOF GLAO mode is GALAGSI WFM. It is ex¬ 
pected to double the ensquared energy within a 0.2” x 0.2” 
sky area at 750 nm with seeing conditions between 0.6” and 


1.1” (Strobele et al. 2012). These partially corrected wave- 
fronts will be sent to the visible integral field spectrograph 
MUSE, located at the Nasmyth B focus of UT4. With the 4 
LGSs pointed 64” off-axis, the MUSE science FoV intended 
to be corrected by this mode is of 1’ across. In this mode, 
the amount of turbulence not detected by the profiler can 
be significant depending on the seeing conditions. It lies be¬ 
tween 2-9 per cent for 2 = 30° for the 14 compiled profiles. 
At zenith, the altitudes can be probed up to 14 km and 
the spatial resolutions in the individual profiles lie between 
252-640 m. 

Fig. 10 shows the performance of the profiler for this 
mode. Here the residuals in CL differ significantly from the 
reference OL ones, with an increase in maximum absolute 
deviations for worse turbulence conditions. Again, we see an 
overexcitation of the first altitude bin in CL, more or less in 
the same range as detected for GRAAL (~10-15 per cent), 
when using the 1240 SAs default mask. This effect in the 
first altitude bin, reflected in the maximum absolute devi¬ 
ation of the residuals, is significantly minimized - down to 
an average of 7 per cent - when the 1044 SAs mask is ap¬ 
plied, in particular for profiles with larger seeings (see the 
dashed thick line in Fig. 10, top-right). However, and as for 
the GRAAL case, applying such mask does not improve the 
result in the second altitude bin for the best seeing profiles, 
this being the reason why the maximum absolute deviation 
of the residuals does not decrease in that regime. This might 
also reflect the need for a more restrictive mask, which would 
decrease the first peak of turbulence and increase the contri¬ 
bution of the second one through the renormalization of the 
profile. In general, the profiler applied to this mode tends to 
identify well all the turbulence peaks above the ground with 
respect to the input profile in the simulations. For the 14 
profiles, the unseen turbulence is retrieved with an average 
and maximum errors of 5 and 12 per cent, respectively. 


4 . 1.3 2-Layer Profile: the F Parameter 

In summary, the profiler succeeds in providing a overall good 
output for the AOF GLAO modes. Both provide unprece¬ 
dented altitude resolution with respect to former systems 
(for instance, the best altitude resolution of GeMS is of 
about 1 km) with extremely good identification of main tur¬ 
bulent layers, and rather fair reproduction of their intensi¬ 
ties. However, it is important to point out that for many ap¬ 
plications robustness is preferred in detriment of spatial res¬ 
olution. Some PSF reconstruction techniques, such as that 
described by Villecroze et al. (2012), for example, only re¬ 
quire a coarse idea of the atmospheric turbulence distribu¬ 
tion (2 or 3 layers). Thus, as an exercise, we can define a 
quantity Fj, as the ratio between the cumulated up to a 
certain altitude h and the total integrated turbulence. Since 
the integration of minimizes the overestimation of the 
intensities in the first two bins, for this test the profiler was 
run with the default mask. 

Running a scan through possible values of h we have 
verified that there is no optimal altitude to compute F/i 
in order to minimize the residuals from the simulations. 
Therefore, ad doc h values of 300 and 500 m are adopted 
for GRAAL and GALACSI WFM, respectively. The re- 
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GALACSI WFM - Open Loop Performance 
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GALACSI WFM - Closed Loop Performance 
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Figure 10. Upper plots: performance of the profiler for GALACSI WFM with the 14 typical Paranal profiles in OL (left) and CL (right). 
Solid thin lines: median absolute deviations; dashed thin lines: maximum absolute deviations; dots: 0 = 0 noise filtering; stars: 0 = 1 
noise filtering; open circles: 0 = 2 noise filtering. Thick lines in the top-right plot refer to results obtained using the mask of only 1044 
SAs and o = 1 noise filtering. Lower plots: two examples of outputs from the profiler using CL simulations: of a good profile with 0.6” 
seeing (left) and of a bad profile with 1.2” seeing (right), both with only ~3 per cent of input unseen turbulence. All output profiles have 
been corrected by the unseen turbulence using as a reference the input value of the simulations, r*Q™. Noise filtering of a = 1 applied to 
the output profiles in the lower plots. 


suiting range of Faoo (GRAAL) is ~0.06-0.86 whereas of 
Fsoo (GALACSI WFM) is 0.15-0.88. The correspondence 
between input (from simulations) and output (from the pro¬ 
filer) F/i parameters is good, with a typical rms error of 
~0.07 for both modes. 

It was interesting to know if the parameter Vh, as de¬ 
fined in the previous paragraph, could be related to the total 
amount of the turbulence corrected by the DSM. For that 
purpose we computed the Zernike coefficients from the POL 
slopes and also from the DSM commands (via slopes space), 
using a proper slopes-to-Zernikes projection matrix with 
150 modes. Possible measurement errors were subtracted 
using the auto-correlation method by Fusco et al. (2004). 
Neglecting piston and TT components, we integrated the 
error-corrected variances of all remaining Zernike modes for 
both POL and DSM-only components. The ratio of the DSM 
variance to the POL variance (called hereafter R) was then 
computed for each profile simulation, resulting in a range of 
0.47-0.87 for GRAAL and of 0.71-0.91 for GALACSI WFM. 
It turns out that there is indeed a correlation of F^ with 
R for both modes, but rather steep, not providing any in¬ 


formation about the altitude until which the turbulence is 
corrected in the GLAO modes. 

An alternative approach consisted in comparing R with 
F = F^/, with h' = min(/i,/iniax), where h is the already 
defined equivalent altitude. Fig. 11 (top) shows how the in¬ 
put (‘sim’) and output (‘cn2’) F correlate for the two AOF 
GLAO modes, in CL regime. Notice that for GRAAL only 
two profiles have h lying inside the probed range of the pro¬ 
filer; so for all the others we integrate the whole profile ob¬ 
tained - properly corrected by the unseen turbulence - to 
obtain F"^. Despite the larger scattering in the results and 
the smaller range probed by this new parameter (0.75-0.92), 
input and output F are fairly well correlated (null-hypothesis 
of correlation with a probability of <0.2 per cent). The pa¬ 
rameter F obtained from the profiler tends to be slightly 
underestimated, with a robust fit giving an angular coeffi¬ 
cient of 0.86 and 0.64 for GRAAL and GALACSI WFM, 


® Given the discretization of the simulated profiles used here this 
ends up being equivalent to integrating the turbulence up to h. 













Tests of a Cf Profiler for the AOF 13 


1.0 

0.9 

^ 0.8 

cQ* 

5, 0.7 

IUh 

0.6 
0.5 
0.4 

0.4 0.5 0.6 0.7 0.8 0.9 1.0 


ft 

/ 






+ GRAAL 
X GAL WFM 
- 1:1 


r (sim) 


1.0 
0.9 
0.8 
lu 0.7 
0.6 
0.5 
0.4 






X + GRAAL (Sim) 

- GRAAL (cn2) 

X GAL WFM (Sim) 
X GAL WFM (cn2) 
- 1:1 


0.4 0.5 0.6 0.7 0.8 0.9 1.0 

^ ~ ^DSM ^ "^POL 


Figure 11. Top: the relation of input T from simulations versus 
output r from the profiler, for the AOF GLAO modes. Bottom: 
r versus the amount of turbulence corrected by the DSM, R. The 
scales of the symbols reflect the amount of turbulence contained 
in the first simulated layer (at 30 m), discretized in 20 per cent 
bins. Both GLAO modes discard the null-correlation hypothesis 
with a probability of less than 1 per cent. 


paring input and output profiles, in particular a criterion 
closer to quantities that matter on the real system (amount 
of turbulence corrected by the AO). 


4.2 LTAO Mode 

GALACSI NFM is the LTAO mode for AOF. This mode 
aims to optimize the correction for the center of the FoV 
- of 20” diameter - and to provide a high correction. This 
translates into a Strehl larger than 5 per cent - with a goal 
of 10 per cent - at 650 nm for a 0.6” seeing and using a NIR 
TT source of at least 15 mag in J-H bands - Strobele et al. 
( 2012 ). 

The main utility of the profiler for this mode is the abil¬ 
ity to locate the most important turbulent layers, useful for 
optimizing its performance. The reconstruction process for 
GALACSI NFM uses a virtual DM approach, described in 
Le Louarn & Hubin (2004). The interaction matrix is cal¬ 
culated using several virtual deformable mirrors conjugated 
to different heights (2 virtual DMs were used in these simu¬ 
lations, conjugated to 0 and 3000 m). These DMs have the 
same number and geometry of actuators as the deformable 
secondary. Using a propagation operator, the cone effect is 
taken into account. Once a synthetic interaction matrix is 
created with all the DMs (virtual and real), it is inverted 
using a MAP algorithm. After that, all actuator commands 
from all the DMs are projected onto the DSM. This opera¬ 
tion is done by simply adding - actuator by actuator - all 
the commands. 

Fig. 12 shows the spatial resolution provided by the HR 
profile, repeating the analysis conducted for other modes. 
Up to 14 km (the maximum height of the simulated profiles) 
about 10 layers can be sampled, and at least 30 per cent of 
the turbulence - considering the 14 typical profiles - will be 
concentrated in the first output bin. The figure shows that 
there is a significant loss of the profiler performance when 
using the default mask; the minimum turbulence detected 
in the first bin, for all profiles, stays always above 81 per 
cent. A noticeable improvement in the residuals is verified 
when using the 1044 SAs mask: the median and maximum 
absolute deviations of the residuals drop from an average of 
~2.8 and 20 to ~1 and 6 per cent, respectively. Even with 
the final configuration for the LTAO mask not yet defined, 
awaiting further improvements on the loop correction, this 
test has shown that by properly masking out problematic 
SAs one can improve significantly the estimation of the 
profile. 


respectively. Assuming a 1:1 relation, the residual rms for 
both AOF modes lies between 0.04-0.06. 

The parameter F correlates very well with R (bottom 
of Fig. 11). An exception occurs for the GRAAL mode with 
profiles that have a very low content of turbulence at the 
ground (profiles ^10 and #13, both with less than 20 per 
cent in the first simulated altitude layer at 30 m). Discard¬ 
ing these two profiles for GRAAL, all subsets shown in this 
specific plot have a null-hypothesis of correlation between F 
and R below 0.1 per cent. Moreover, assuming a 1:1 relation, 
the rms of the residuals lies between 0.03-0.06. In summary, 
R seems to be a good indicator of the correction performed 
by the DSM up to the average altitude h, as long as R> 0.6. 
Also, the parameter F provides a better criterion for com- 


4.3 Final Remarks 

It has been verified that for the 3 AOF modes - GRAAL, 
GALACSI WFM & GALACSI NFM - it is possible to char¬ 
acterize relatively well the turbulence distribution with the 
chosen profiler once a correlation mask is properly de¬ 
fined. The importance of this mask seems to increase with 
the level of correction expected from the mode, being cru¬ 
cial for LTAO. The major peaks of the input distribution 
are identified in all modes, and quantitatively, at least for 
the 1044 SAs mask chosen here, the median absolute devia¬ 
tions of the profile residuals lie within ~2.5 per cent. There 
is still some intensity cross-talk between the very first bins. 
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Figure 12. Upper plots: performance of the profiler for GALACSI NFM with the 14 typical Paranal profiles in OL (left) and CL (right). 
Solid lines: median absolute deviations; dashed lines: maximum absolute deviations; dots: a =0 noise filtering; stars: a =1 noise filtering; 
open circles: a =2 noise filtering. Thick lines in the top-right plot refer to results obtained using the mask of only 1044 SAs and a =1 
noise filtering. Lower plots: two examples of outputs from the profiler using CL simulations: of a good profile with 0.6” seeing (left) and 
of a bad profile with 1.2” seeing (right). Noise filtering of a = 1 applied to the output profiles in the lower plots. 


specially in the best resolved case - GRAAL - which has 
to be further addressed if such a fine resolution is actually 
needed; this cross-talk is the major cause of the large max¬ 
imum absolute deviations of the residuals observed for this 
mode (up to 20 per cent). 


5 SUMMARY AND CONCLUSIONS 

We have validated the so-called wind-profiler algorithm un¬ 
der the AOF configuration. This validation was carried out 
using end-to-end Octopus simulations. 

In a first instance, OL simulations allowed us to check 
the profiler limitations and to tackle its sensitivity issues in 
aspects which are intrinsic to the profiler, such as the effect 
of uncertainties in the input parameters (the outer scale and 
Na layer altitude), sampling effects and response to different 
types of turbulence distribution. 

The detected noise level in bins with simulated zero tur¬ 
bulence is consistent with 0.5 per cent error. We have, how¬ 
ever, found some leakage from the profile output intensities 
whenever the simulated turbulent layer is strong and falls 
exactly on an output bin altitude. This amount of leakage. 


which can reach up to 20 per cent for a single input turbulent 
layer, can be regulated to some extent by adopting a noise 
filter. When spatial resolution is not an issue, the integra¬ 
tion of neighbouring output bins can decrease the residuals 
with respect to the input turbulence. Alternatively, it is also 
possible to retain only the LR output profile instead of the 
combined one. 

The absolute calibration of the output profile depends 
strongly on an assumption for the turbulence outer scale, a 
parameter which cannot be directly retrieved with the pre¬ 
sented algorithm. A way to circumvent this problem is to 
get an independent estimate of turbulence parameters and 
to use them in the assumed atmospheric model for the pro¬ 
file. Alternatively, it is possible to opt for the method pre¬ 
sented in Valenzuela et al. (2014), which makes an estimate 
of the total Lq by fitting the slopes autocovariances and 
highlights the need for a more detailed study of the outer 
scale altitude distribution in order to properly calibrate the 
profile. Bearing these calibration issues in mind and con¬ 
sidering that the most important information for the AOF 
performance improvement is contained in the relative 
distribution, our work was focused on testing only normal¬ 
ized profiles. Nevertheless, Fried parameter estimates given 
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by the profiler and by the WFSs slopes variances could still 
be useful for the correction of the unseen turbulence, regard¬ 
less of the assumed outer scale. 

Dataset length and undersampling can have an impact 
on the normalized estimation. This study recommends 
for the AOF a minimum acquisition of 15k data frames to 
feed the SPARTA profiler routine, which is then submitted 
to an undersampling of r = 10. Preliminary tests running 
in a parallelized Matlab environment have shown that the 
profile computation of such amount of data is viable in about 
15 seconds, and therefore suitable for real time estimations. 

In CL, a reliable estimation largely depends on the 
knowledge about the quality of correction achieved all over 
the pupil. This happens in particular for GALACSI NFM, 
and, to a lesser extent, to GALACSI WFM. The design of 
realistic correlation masks - that in the future should take 
into account features such as obscuring spiders and possible 
defective actuators - are a crucial step in order to avoid the 
introduction of spurious correlation signals in the profiler. 
The masking procedure has already been used in previous 
works with this profiler (e.g. for GeMS), but might need a 
special attention here because of the different DM nature. 

Using a symmetrical ad hoc mask of 1044 SAs to anal¬ 
yse simulations with the 14 typical Paranal turbulence pro¬ 
files, the profiler has provided good estimates for all the 
AOF modes. The GLAO modes are able to provide an un¬ 
precedented spatial resolution in the determination. It is 
recommended that unseen turbulence correction of the de¬ 
rived profile is carried out only for the GRAAL mode. 
For GALAGSI modes this correction is either unnecessary 
(NFM) or negligible (WFM), in particular for the latter if 
we consider the uncertainties involved versus the amount 
of unseen turbulence considered. If a fine resolution for the 
profile is not required in the GLAO modes, it is possible 
to define a parameter F, integrating the output profile 
up to the altitude h, determined from the same profile. The 
parameter F has then a direct relation with the amount of 
turbulence corrected by the DSM. The profiler used with the 
LTAO mode succeeds in providing the locations of the main 
turbulence peaks with a spatial resolution within ~ 1.0-1.5 
km, enabling one to fine tune the loop reconstruction process 
in real time. 
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